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GENERALIZED TREATMENT OF PLANE ELECTROMAGNETIC WAVES 
PASSING THROUGH AN ISOTROPIC INHOMOGENEOUS PLASMA SLAB 
AT ARBITRARY ANGLES OF INCIDENCE 
By Calvin T. Swift and John S. Evans 


SUMMARY 


The differential equations describing the passage of plane electromagnetic 
waves through a plane -parallel inhomogeneous plasma slab at an arbitrary angle of 
incidence are solved by numerical integration between appropriate boundary con- 
ditions. The methods used are described and a typical solution is discussed. 


INTRODUCTION 


Analytical solutions to the wave equation in an inhomogeneous plasma can be 
derived for simple electron density distributions (refs. 1, 2, and 3) j "but such 
solutions are often inadequate for application to a given practical problem. One 
approach (ref. 4) considers the plasma as a multilayered series of homogeneous 
slabs. This technique has been derived for convenient application to waves nor- 
mal to the plasma slab; however, extension to include oblique incidence seems 
questionable because derivatives of the plasma properties which generally appear 
in the propagation equations for inhomogeneous media do not explicitly appear in 
the multilayered-slab approach. The W.K.B. approximation (refs. 5 and 6) gives 
good solutions for arbitrary electron density distributions so long as the con- 
ditions for its validity are met. These conditions are not met if plasma prop- 
erties vary appreciably within a wavelength; and since such plasmas are fre- 
quently of interest, another method for obtaining solutions is needed. 

The purpose of the work reported herein is to reduce the propagation equa- 
tions to a form suitable for numerical integration on an electronic computing 
machine. The physical model postulates arbitrarily polarized, plane, monochro- 
matic electromagnetic waves impinging at an arbitrary angle of incidence on a 
slab of plasma bounded by infinite parallel planes, which separate the plasma 
region from vacuum. The plasma properties are assumed to vary only in the direc- 
tion of the normal to the planes. In general, discontinuities in the plasma 
properties and their first derivatives are permitted only at the boundaries; how- 
ever, at normal incidence, the continuity requirement on the first derivative can 
be removed. 



SYMBOLS 


In all derivations mks units are used, but in tables and graphs, linear 
dimensions are shown in centimeter units. 


c 

E 

F y 

f 

% 

if 

-» 

k 

M,N 

N e 

n 


T 1^ T 2 


Sc 


Uu 


speed of light, 3 X 10® meters per sec 
electric field intensity, volts/meter 

z-dependent part of y-component of electric field intensity in plasma, 
volts/meter 

frequency, sec*" 1 

z-dependent part of y-component of magnetic field intensity in plasma, 
ampere -turn s/me ter 

magnetic field intensity, ampere-turns/meter 
propagation constant, o>/c, meters" 1 

propagation vector, k x ii x + k z u z = k q sin 0ii x 4- k q cos 0u z , meters* 1 

arbitrary constants determined by solution of wave equation, volts/meter 
number of electrons per cubic meter 


index of refraction. 


W 


voltage reflection coefficients of perpendicular and parallel components 
of wave, respectively 

position vector, xa x + ySy + zii z , meters 

voltage transmission coefficients of perpendicular and parallel compo- 
nents of wave, respectively 

unit vector in x-direction 
unit vector in y-direction 


u 9 


unit vector in z-direction 


ii ■■ ii 



V = 1 - 


1 



w = 


v/tup 1 



x,y,z 


X 

Z 

z 0 

Z 1 


7 

Sr 

*t 

€ 

G o 

0 

A 

^0 

v 

I 

0 


Cartesian coordinates 

x-dependent part of solution of wave equation 
z-dependent part of solution of wave equation 
thickness of plasma slab, meters 

ramp depth for trapezoidal electron density, meters 

amplitude factor defined by equation (110), radians 
phase angle of reflected wave, radians 

phase angle of transmitted wave, radians 
permittivity, farads/meter 

permittivity of free space, 8.854 X 10“^ farads/meter 

angle of incidence, deg 
angle of transmission, deg 

wavelength, meters 

permeability of free space = 4jt X 10 “7 henry/meter 

mean collision frequency of electrons with neutral particles, 
collisions per second 

phase difference between parallel and perpendicular components of 
incident wave, which accounts for elliptical polarization, radians 

angle between incident electric vector and y-axis, deg 
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i|r phase factor defined by equation (ill), radians 

cjd frequency of propagating wave (radians per second) 

CDp plasma frequency (radians per second), 2ir X 8.97 X \j^e 


j j magnitude of a complex number 

T differentiation with respect to z 

Subscripts: 

E components of X and Z which are solutions to the electric field 

intensity 

fs wave in free space (includes incident and reflected components) 

H components of X and Z which are solutions to the magnetic field 

intensity 

i incident wave 

p wave transmitted into plasma 

r reflected wave 

t wave which has transmitted through plasma slab 

x x-component of wave 

y y-component of wave 

z z -component of wave 

Tfs tangential component of wave in free space (including incident and 

reflected waves) 

Tp tangential component of wave which is transmitted into plasma 

0 amplitude of wave or free-space parameter, as appropriate 

00 amplitude of arbitrarily oriented incident wave 

1 component of the wave which is perpendicular to plane of incidence 

2 component of the wave which is parallel to plane of incidence 

Note that the subscripts are used in combination. For instance, E r0 y is 
the amplitude of the y-component of the reflected electric vector. 
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THEORY 


Electromagnetic Waves in an Inhomogeneous Plasma 

— > — > , . 
The wave equations for E and H in inhomogeneous media are (refs. 1 and 7) 


V^E + kQ^n^E = 



( 1 ) 


+ kQ 2 n 2 H = 



( 2 ) 


where the index of refraction n is a function of position. 


Let an electromagnetic wave he incident upon an inhomogeneous plasma slab 
as shown in figure 1. The plane of incidence of the wave is defined as the plane 
which contains u z and T sq . The 
electric vector is, by definition, 
orthogonal to kQ; however, there is 


no imposed restriction on the orienta- 
tion of I? with respect to the plane 
of incidence. Although a ray path is 
shown in the diagram for visual pur- 
poses, the method presented does not 
depend upon explicit knowledge of the 
ray path. It is of interest to note 
(see appendix) that the angle of the 
emerging ray is identical to that of 
the incident wave, even for the inho- 
mogeneous case discussed in this anal- 
ysis. If the variation of the index 
of refraction is restricted to one 
dimension so that 

n = n(z) (3) 



Figure 1.- Geometry of plane-wave inter- 
action with an inhomogeneous plasma slab. 
N e = N e (z)> v = v(z). 


and. if the wave in free space is independent of y , then the solutions to equa- 
tions (l) and (2) are independent of y. (See appendix.) If the condition pos- 
tulated by equation (3) is used and all derivatives with respect to y are 
ignored, equations (l) and (2) become 


-2-» 
d E 


n 2-> 
d E 


^ + ^ + k^ n 2( z )E = --2. 

x 2 ^ 2 ^ 
ox oz 


a 

E dn 2 U) 


E an 2 U> 

chc 

Sz 

u - — 

x a z 

% Sz 




(M 
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~ 2 — > ^2 — > » 

ll. + IS + k 0 2 n 2 (z)H = 
ax 2 8z 2 n2 (z) 


1 8n ( z ) / dH x 

dz \8z 



1 5n 2 (z) 


n 2 ( z ) Sz 8z 


— (5) 


Solutions to equation (4) for arbitrary orientation of E can be obtained 
by solving three differential equations in which E x , Ey, and E z are dependent 

variables. The solutions of equation (5) can be found from the solutions of 
equation (4) by means of Maxwell’s equations. However, Stratton (ref. 8, p. 492) 
has shown that the wave in the plasma can be uniquely defined by two rather than 
three differential equations in which the dependent variables are E]_ = Ey and 

E 2 = \|e x 2 + E z 2 and, since E 2 is related to Hy through Maxwell’s equations, 
the two differential equations for which solutions are required are: 


For the Ej_ case: 


For the E 2 case: 


th. + th. + , 0 


( 6 ) 


8x' 


8z‘ 


th 

ax 2 


+ + ko 2 n 2 (z)H y = — i— h\ z ) 5k 

8z 2 7 a2(z) Sz 8z 


(7) 


If a product solution to equations (6) and ( 7 ) is assumed, it follows that 


Ey(x^z) = Z E (z)X E (x) (8) 

and 

Hy(x,z) = Z H (z)X H (x) ( 9 ) 

The substitution of equations (8) and ( 9 ) into equations (6) and ( 7 ) results in 
a solution for X E and Xg. Solving for the X's yields for Ey and Hy 


Ey(x,z) = Fy(z)e ± ik x x ( 10 ) 

Hy ( x j z ) = Gy(z)e ±i ' li x x (11) 
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But from the arguments of the appendix and from figure 1, k x = kQ sin 0- If 
the condition is made that the wave propagates in the positive x-direction, the 
negative sign on the exponent can he dropped. The positive sign on k x x implies 

also a time factor of e"^ 0 ^. Equations (lO) and (ll) now become 

Ey(x,z) = F y (z)e lk O XSin9 (12) 

H y (x } z) = Gy(z)e lk ° XSin0 (13) 

Substitution of equations (12) and (13) into equations (6) and (7) further 
reduces the wave equations to 


Fy [ Z '~ + ko 2 [p 2 (z) - sin 2 e]Fy(z) = 0 (l4) 

dz^ 


a 2 Gy( z ) 

dz 2 


dn 2 (z) dGy(z) 
n 2 (z)dz ^- z 


ko 2 [n 2 (z) 


sin 2 e]Gy(z) = 0 


(15) 


The solution to the wave equation in free space is obtained by setting 
n 2 = 1 and the right-hand side of equations (l) and (2) equal to zero. If the 
propagation vector kQ lies in the xz-plane, then the incident waves will he of 

the form 


E. = e* e ik 0i‘ r 

^iy ^iOy e 


%y = HiOy* ill0i ' r 


(16) 


From inspection of figure 1, it can readily he determined that 


kQ^ • ~r = kQ^sin 8u x + cos 0u z ) • ^xu^- + yu^- + zu z ) 


\ " “/ \ 
= kQX sin 0 + Fqz cos 0 


(17) 
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Likewise, the reflected waves have the form 


with 




E ry = E r0y e 


ikor* r 


H ry = H r0 ye ik Or ' ] 


(18) 


— ^ ^ 

k Q r * r = kgx sin 0 - kgz cos 0 


and the transmitted waves will be 


Et y = E t 0 y e ik Ot*r 




H t y = 




(19) 


( 20 ) 


with 


kQ^. • r = kQX sin 0 + kQZ cos 0 


( 21 ) 


Boundary Conditions and Reduction of Wave Equations 
to Computational Form 

Case 1: E perpendicular to the plan e of in c idence .- When E is perpendic- 

ular to the plane of incidence, the boundary conditions at the surfaces z = 0 
and z = Zq require that the tangential components of if and if be continuous 

across each surface. At the boundary z = 0, 


E TfS ~ E Tp 

(22) 

frrfs = H T p 

(23) 


Equation (22) requires that 

X fooyUy + Er0 ySy)e 11S 0 XSlne - S, X SyFy( 0 )e ik O XSin0 (24) 
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or 


No v, since 


equation (2j) becomes 


E i0 y + E rOy - F y(°) 


— » q — > —* 

H = — — V X E 
i(X4 0 


(25) 


(26) 


uyk Eiy + ^ ^ ' %e lls 0 x5ln %’ (0) (27) 

or 

ik© cos 9 ^-®i0y ~ E r0y^ = E y (®) ( 28 ) 

Similarly, the boundary conditions at z = Zq lead to the equations 

Fy(z 0 ) = Et 0 ye ik ° Z0COSe (29) 

and 

Fy’(z 0 ) = iko cos 0E tO ye ik O z O cose (30) 

Thus, the problem is reduced to solving the differential equation (l^): 


y ( ^ + ko 2 [n2(z) - sin 2 e]Fy(z) = 0 

dz 2 

subject to the four boundary conditions (eqs. (25), (28), (29), and (50)). The 
wave equation and the boundary conditions are composed of complex quantities. It 
is therefore necessary to separate real and imaginary parts prior to machine pro- 
graming. If the boundary conditions are divided by E^Qy a nd if the following 

definitions are expressed: 


E r0y 

E i0y 


|Rl|e 15 rl 


(51) 
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then the boundary conditions (eqs. ( 25 ), ( 28 ) ^ ( 29 ), and (30)) reduce to 


%0y 


LKq cos 0 1 - Rj_ 


•iSrx'l - V<£> 


F y (z 0 ) _ I I itkozocose+^tl) 


= Tl e 


F y'( z o) 


ikg cos 0 Tj_ 


.[kQZocose+Stl 


If the boundary conditions are divided by 


i( k OZO cos0+5 tl) 


and if further notations are introduced, such that 


- i ( k O z O cos9+5 tl) 


- _ 1 R 1 I -i(koz o cos0+5 tl ) 


R-i = e 

u I _ 1 


i 




the boundary conditions (eqs. ( 33 ) to ( 36 )) now become 


Ii - V Srl - ^ (0) n 

EiOy 

(39) 

lk 0 cos - R l el6rl ) - E 10y 

(40) 

^(zo) 1 ! _ ± 

(41) 

E iOy 


F y'( z o)ii 

= iEq cos 0 

E iOy 

(42) 

If the differential equation (l4) is multiplied by Ii^Ej_oy> 

the normal! za- 

tion of the equations is complete. If the following definitions 

are expressed: 

Fy(z)l]_ , , . , . 

— = r(z) + is(z) 

EiOy 

(43) 

n 2 (z) = V(z) + iW(z) 

(44) 

V*(z) = V(z) - sin 2 0 

(45) 


the differential equation (l4) becomes 


d + 1S ^ + kQ 2 (V* + iW)(r + is) = 0 (46) 

dz 2 

Equation (46) can be separated into real and imaginary parts , and they must 
independently vanish; therefore. 
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(^ 7 ) 


+ ko 2 (V*r - Ws) = 0 

dz^ 

+ kQ 2 (Wr + V*s) = 0 

dz 2 

Also, n^(z) is Interpreted as the square of the complex index of refraction 
with 


V = 



W 



From the definition (43) > the boundary conditions at Zq (eqs. (hi) 
and (42)) become 


r ( z o) 
s ( z 0) 



r'(z 0 ) = 0 
s’(zq) = ko cos 





m 


(^ 9 ) 


The parameters Ip and Rp can he solved from the boundary conditions at z = 0, 


Ip + Rpe 1 ^ 1-1 = r(o) + is(o) 


(50) 


ikQ cos 9^Ip - Rpe'*’^ 1 ’ 1 ^ = r'(o) + is' 


( 0 ) 
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( 51 ) 


The solutions to these quantities are: 


and 


x l = | 


'( 0 ) + 


s‘(0) 


k q cos 0 


+ ii 
2 


s (o) - 


( 0 ) 


k q cos 0 


= |(Ai + iBi) 


( 52 ) 


Ri = 


r (0) 


s‘(0) 


kQ cos 0 


+ ii 


b(0) + 


r«(0) 


kQ cos 0 


-iS r i 


§(Cq + iD 1 )e“ l8rl 


(53) 


The values r(o) , r’(o), s(0), and s'(o) are found from the solution of 

the wave equation ( 47) . 

From equation (37); the magnitude of Ij is 


!l| = 


Til 


(54) 


Therefore, 


I Til = 


\I A l 2 + ®1 2 

And from equation ( 38 ) , the magnitude of R! is 


(55) 


R 1 


I R 1 

|Ti 


(56) 
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Therefore, the reflection coefficient is 


l R il 


C ± 2 + Dj 2 
\j A X 2 + Bq 2 


(57) 


The quantity 5 r q is the phase of the reflected E vector and S^l is the 

— ^ 

phase of the transmitted E vector. The phase factors are computed by first 
making equation (53) equal to the definition of equation (38)* Thus, 




l R il 

Ti| 


^-i^kozocose+Stl) 


But, from equation (56) 


(58) 


Ri 


Ri 


1 I = 


T 1 


- 


2 + Dq 2 


Therefore , 


C 1 + iD l _ -i(kozo cos0+s tl" 5 rl) 
\[ Cl 2 + Dl 2 


and 


or 


tanfkQZQ cos 9 + 8 t i - S r q) = - ~ 


/ D P 

EqZq cos 0 + 6-^1 - 8 r i = arc tan - — 

\ c i t 


Similarly, if (52) and (57) are equated, 



i(kozocos0+5tl) 


I Tl| 


(60) 


(61) 


(62) 
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But, from equations (5^) and (55) 


Therefore, 



1 

2 


\j Aq 2 + Bj 2 


A 1 + -i^koZQCOsS+S-ti) 

\) A l 2 + B x 2 


and 


kQZQ cos 0 + 5-^1 = arc tan - 



(65) 


(64) 


Equations (6l) and (64) are equations with which to solve for the two unknowns 
S r l and &tl» The solutions for the phase factors are: 


B 1 

B-tl = arc tan[ J - kQZQ cos 0 

A 1 


(65) 


/ M / 

^rl = arc tan I - — - arc tan - 


(66) 


Case 2: E parallel to the plane of incidence .- As in case 1, the boundary 

conditions at the surfaces z =_0 and z = z 0 require continuity of the tan- 
gential components of E and H. 

At z = 0, equation ( 23 ) establishes the requirement that 

u z X ^H^oy u y + B r0y u y) = u z x Gy(o)uy ( 67 ) 

or 

B i0y + B r0y = ^y(^) (67 a ) 
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And, since 


— ■ > T — > — > 

E = — V X H 

iCD£ 


equation (22) requires that 


u z x (- -7-=— V x H ± - 7-^— V x H r ) = 


icoe 0 


xcoe 0 


-u'z X V X Gy(o)uy ikQxsin0 
io>e(o) 


or 


ikQ cos 0(H iOy - H r0y 


G v ’(0) 


n 2 (o) y 

Similarly, the houndaiy conditions at z = z q require that 


^ , . TT ikoz n cos0 

G y( z 0) _ HtOy e 


and 


G y’( z 0) = n2 ( z 0) ik 0 cos 0HtOy elk ° ZoCOS0 


As in case 1, the problem is reduced to solving equation (15) 


d 2 Gy(z) p dn 2 (z) dGy(z) 2 To/ \ o"l / \ 

-tr - -sH- + **> - sln2e >^> - 


0 


dz £ 


subject to the boundary conditions (67a), (69a), (70), and (71)- If the 
conditions are normalized as in case 1, they reduce to 


I 2 + R 2 e l5r2 = i 2 = t(o) + iu(o) 

H iOy 


ikQn 2 (0)cos 0 /l 2 - R 2 e i5 r2^ = JL - l i 2 = t'(o) + iu'(o) 
' 7 H iOy 
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( 68 ) 

(69) 

(69a) 

(70) 

(71) 

boundary 

(72) 

(73) 


^ 12 = t(ZQ) + lu(zo) = 1 
H iOy 

(7*0 

Gy* ( z o) 

g- 0y 12 - t (zo) + iu( Z0 ) 


= ikQn^(zo)cos 9 


= ikoV(ZQ)cos 0 - kQW(zo)cos 9 

( 75 ) 

n 2 ( 0 ) = V(o) + iW(O) 

( 76 ) 

n 2 (z 0 ) = V(z 0 ) + iW(z 0 ) 

( 77 ) 

I* 0 * . | R 2 |e iS r2 
H iO y 

( 78 ) 

•“iOy 

(79) 

-i(koz o c°s0+8t2) 
12 = 1 

l*a| 

(80) 

* _ t R 2 1 - 1 ( k 0 z 0 cose+5 t2) 

ilp ' G 

|t 2 | 

(81) 

J CVJ 
H Eh 

II 

OJ 

H 

(82) 

|r 2 | = l%l 1 
l T 2 1 

(83) 
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12 = t(z) + iu(z) 


Gy( z ) 

H iOy 


(84) 


If the differential equation (15) is multiplied by l^f^lOyf the following equa- 
tion results: 


i-Ytlz) * iu( z fl . &(f) - y < z .fl /fghl M£) 
dz 2 |v 2 (z) + w 2( z )] \L dz dz 


dW(z) du( z )~ | + i [~ dW(z) dt(z) 
dz dz dz dz J I dz dz 


dY ( z ) du( z ) 
dz dz 


I 


)] + ijw 


+ kQ < V*(z)t(z) - W(z)u(z) + i W(z)t(z) + V*(z)u(z) > = 0 


4} 


(85) 


As before , the real and imaginary parts must independently vanish, whereby two 
simultaneous differential equations result: 


d 2 t ( z ) 

V(z) 


dV(z) dt(z) dW(z) du(z)~j 

W(z) Id 1 

dz 2 jy 2 ( 

>) + W 2 (z)J ^ 

_dz dz dz dz J 

r 2 (z) + W 2 (z)j *- 

! dV(z) du(z] 
dz dz 


V*( z 

-)t(z) - w(z)u(z)j = 0 



dW(z) dt(z) 
dz dz 


(86a) 


d^u(z) 

V(z) 

IdW(z) dt(z) + dV(z) du(z) 

+ - 

W(z) 

IdV(z) dt(z) 

dz 2 

V 2 (z) + W 2 (z) 
_ _ 

L dz dz dz dz _ 

[v 2 ( z) + w 2 (z) 

[_dz dz 


- Mil Mil + kQ, 2 W(z)t(z) + V*(z)u(z) = 0 (86b) 

dz dz 


By separating real and imaginary parts in equations (74) and (75) the 
boundary conditions for t and u at zq are: 
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(87) 


^ZO) 

u ( z o) 



^(zq) 

u '(z 0 ) 


= -kQW(zo)cos e 
= 15qV(zo)cos 0 


( 88 ) 


The boundary conditions, and the solutions of u and t at z = 0, 
determine the reflection and transmission terms. 

From equations (72) and (73) 


and 


I 2 = 


||/t(0) + 


-w(o)t'(o) + V(o)u , (Q) 
k 0 COS e[w 2 (o) + v 2 (o)lj 


+ i C(o) + -v(o)t'(o> w(o)u' (o) 

| k q cos e[w 2 (o) + v 2 (o[| 

= |(a 2 + iB 2 ) 


R 2 e i6r2 = | f ft(0) + W ( 0)t»(0) - v( 0 ) u'(0) 

1 | kQ cos e|w 2 (o) + v 2 (o)] 

+ i t(o) + ^0)^(0) w(o) _u'(o) 

\ ko cos e[w 2 (o) + v 2 (o)]j 

- |(°2. + ^ 2 ) 


(89) 


(90) 
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The transmission and reflection coefficients are: 


l*2| 


1 _ 2 
1 12 | | a 2 2 + B 2 2 



| r 2 

1 A 

c 2 2 + d 2 2 

\M 

II 

a 2 2 + B 2 2 


( 91 ) 




( 92 ) 


Equations ( 89 ) and ( 90 ) (with definitions ( 80 ) and (8l)) are equivalent in form 
to the corresponding equations for E perpendicular to the plane of incidence; 
therefore, the phase terms of the parallel component will be 


5 t2 


= arc tan - 



- JsqZq cos 9 


( 93 ) 


6 r 2 


= arc tan - 



arc tan - 



( 94 ) 


Arbitrary Orientation of Electric Vector With Respect to 
the Plane of Incidence 

If the electric vector is oriented at an arbitrary angle with respect to the 
plane of incidence, it may be resolved into vector components parallel to and 
perpendicular to the plane of incidence. The purpose of this section is to see 
how the components recombine after the wave interacts with the plasma. 

Case 1: Transmission coefficient for E at arbitrary orientation to plane 

of incidence .- From figure 1, let be the angle between the electric vector 
and the y-axis. If the amplitude of the electric vector is Eqq> then 


E iOy = E 00 cos 0 

EiOx = "EOO sln 0 cos 9/ 

E iOz = %) sln 0 sin 9 ^ 


t 


( 95 ) 
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and since (eq. 26 ) 


H = -r±— V X E 

ioho 


H iOy 


_ 1 ^Eiz j e -i(koxsin0+k0zcos6) 

icqiQ y 6z 8x J 

[5- E 00 sin 0 
^0 


the transmitted waves are: 


E ty - ITJEoo =o« 0e i ( k O xsin0+k O zcos0+s t 1 ) 
= ®tl 


H ty = “\/^ 1T 2 |E 00 sin 
' ' 4 q 


^i(koxsin0+kozcos0+6-( ; 2) 


But 


-f' 


E t2 _ UEtx^ + E tz^ 


where E^ x and E-^ z can be determined from the relationship (eq. (68)) 


-i — > — > 

E — V X H 

iCD€n 


Performing this operation yields 

St2 - (-cos e? x + ain e^)| am )Z ie 1 ( k o xsine+koz<=ose+&t2 ) 

and, from equation (99) 

E t2 - |T 2 |Eoo sin j fc i(‘b”in«+M“»8 + 6 t2 ) 

[U o 


= -H 


ty 


(96) 

(97) 

(98) 

(99) 


( 100 ) 

( 101 ) 
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The tr ansmi tted electric vector is the vector sum of its components, or 


E t = \I E tl 2 + E t2 2 


(102) 


Substituting equation (97) and equation (101) into equation (102) and dividing 
by Eqq gives 

\!/2 


E. 


— — = ( Ti 2 cos20e 

% \ 


2i5tl + T 2 2 sin20e 215t2 J e i(k> xsin 9 +k ozcose) 


1/2 


Ti 2 cos 20 + T 2 2 sin 2 0e 2i ^ t2 " 5tl ) e i(k o xsin0 + k oZ cose + 5 t i) 


(103) 


Let 


2 2A m2 • 2A 2i (^b2“^tl) 
Tq^cos^ + T 2 ^sin c -0e 


1/2 


|^|ti 2 c o s 20 + T 2 2 sin 2 0 cos 2^ 2 - S^q) 

it 


+ i 


T 2 2 sin20 sin 2 (St2 “ 5 tl) 


= ye 


and define 


1/2 


(104) 


= Tq 2 cos 2 0 + T 2 2 sin 2 0 cos 2(&t2 - &tl 


Then 


= T 2 2 sin 2 0 sin 2^&t 2 - &ti 
Aj + iBj = 7 2 e 2 ^ 

= 72 cos 2 i|f 

B^ = 72 sin 2t 


(105) 


(106) 


(107) 

(108) 

(109) 
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A 5 2 + B 5 2 = 


1 / b 3 
ilr = — arc tan — — 

2 \ A 3i 


Therefore, 


- ( a 3 2 + b 3 2 ) e 


i kQXsin0+kQZcos6+5^^+ — 


— arctan — 


i t i - §*- « (a 5 2 + b 5 2 


It may he noted , in general , that a polarization shift is implied hy 
equation (112). 

Case 2: Reflection coefficient for E at arbitrary orientati on to plane 

of incidence.- As in case 1 (eqs. (95) and (9&J) 


E iOy - E 00 cos 


~ E oo sin 0 


But from equations (l8), (19); (51 ) > (78), (95); and ( 96 ), the reflected com- 
ponents are: 


E ry - | Hi | Boo GOS 0 e 


i ^kQXsinG -kQZcos0+5 r i^ 


e if kQXSin0 -kQZCOs0+5 r 2 

Hr y = -v | ^2 | Eqo sin 0 e 

\ ^0 


From the relationship (eq. (68)), 


E = - - V X H 

10*0 
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equation (115) determines the parallel component of the electric vector. When 
equation (68) is used. 


-» / _» A. , . i(k 0 xsine-k 0 zcose+6 r2 ) 

E r2 - ( cos 9 u x + sin 0u z/ 1 R 2 |E00 sin 0 e 


(116) 


and 


2 2 

E r 2 = \jErx + E rz = | R 2 | E 00 sin 0 e 


Combining equations (ll4) and (117), with the definition 


i ^QXsinO -kQZCOs0+S r 2^ 


(117) 


E r = 1/ E rl 2 + E r2 2 


(118) 


and dividing hy Eqq yields the reflected component 


E 00 


r = ( Rj_^cos^0e r "*" + R 2 ^sin 


2i5 r 2^ 1//,2 e i ( k O xsin0-k O zcos Q) 


Rj^cos^ + R2^sin20e 


2i(s r2 -6 rl ) 


1/2 


1 (^0 xsin6 -fcozcos0+5 r q) 


(119) 


The form of equation (119) is similar to that of equation (103). Therefore, 


= (ai , 2 + Bi, 2> ) 


=“00 


1/4 1 


A + B Ifj e 


/ B \ 

|^kQXsin0-lQZcos0+6 r q+ ^ arctanf^jy 


( 120 ) 


I R I = 


E, 


E 


'00 




+ Bb 2 ) 


1/4 


(121) 


where 


A4 = R^ 2 cos 2 0 + R 2 ^sin20 cos 2 - 6 r i 

B4 = R 2 ^sin20 sin 2^5 r2 - &rl) 


( 122 ) 

(123) 
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Arbitrary Polarization of Incident Wave 


The electric vector of an elliptically polarized wave can be resolved into 
two orthogonal linearly polarized components. Suppose the incident wave of fig- 
ure 1 is elliptically polarized so that 


-> _ -> ik^r+ij- 

Eiy - Ej_QyUye 

where £ is the phase difference between the perpendicular and parallel compo- 

i£ 

nents. If E^Qy is replaced by E iOy e > the definitions of the reflection and 
transmission coefficients are altered so that equations ( 31 ) and ( 32 ) become: 


E rOy 

E 10y 



(124) 


Iflle 1 ^) (125) 

E iOy 


The boundary conditions (33) "to (3 6 ) are also modified so that they become: 


X + R^e " -^l 


Fy(0) 

E 10y 


Ine -i £ 


( 126 ) 


ikQ cos 0(li - R]_e i5r l 


V(°) 

E iOy 


In e 


-ii 


(127) 


Fy — = 1 (128) 

E iOy 


— = iko cos 6 ( 129 ) 

E iOy 
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However, the additional factor § can be absorbed by redefining r and s 
(eq. (43)) 


r(z) + is(z) 


Mil I e-l! 

E iOy 


(130) 


If the wave equation is multiplied by — — e"^, then the differential equa- 

E iOy 

tions which involve r(z ), s(z), and the boundary conditions remain unaltered. 

Thus, the appropriate parameters for an elliptically polarized wave are obtained 
by adding the phase factor | to the solutions for 8 r p and 8ti* 


RESULTS 


Test Cases for Machine Program 

The numerical integration of the propagation equations was performed on an 
IBM 7090 electronic data processing system by the Runge-Kutta method with an 
accuracy of 10“7 per integration. 

As a first check of the computer program, analytical and machine calcula- 
tions were compared for the homogeneous collision-free plasma (v jib p = oj. The 

value of — = (hence an index of refraction of l/l.5) was chosen because 

“p fi 

Rossi (ref. 9> pp. 137 an< i 377) has computed convenient parameters for n = 1 / 1 . 5 , 
which aided hand calculations. The plasma slab was arbitrarily chosen to be 
6 free-space wavelengths deep. Figure 2 shows that the agreement was good. 

Machine computations were also carried out for cases with collisions in order 
to observe the behavior of the transmission coefficient when losses are introduced. 
The results of these calculations are in figure 2. For the lossy plasmas, hand- 
calculated transmission coefficients at 0=0 (normal incidence) were in excel- 
lent agreement with the machine results. The interference effects, which appear 
so strongly for v/oip = 0 are caused by multiple internal reflections. As the 

plasma losses become more predominant the waves which undergo multiple reflections 
are damped out, as figure 2 shows. Figure 2 also shows that the calculated trans- 
mission coefficient drops sharply as 0 exceeds the critical angle, which is 
defined as the angle of incidence for which the refraction angle is 90°. For 
angles of incidence greater than the critical angle (0 « 42°) the transmitted 
wave is exponentially damped (ref. 8, p. 498). 

Figure 3 is a plot of phase difference between the reflected y-components of 
E and H as a function of the angle of incidence. Rote that figure 3 represents 
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Figure 2.- Transmission coefficients for a homogeneous plasma slab. 


the phase difference between 




provided l8o° is 


added to every point on the ordinate. An interesting feature of this plot is the 
abrupt change in phase difference at 9 ~ 34°, which is the Brewster , or polar- 
izing, angle for this medium (ref. 8, p. 497). 


Table I lists all of the computed results obtained for the homogeneous plasma 
with — = — . 

“p p 


Plots of reflection and transmission coefficients at normal incidence for 
a plasma slab with constant collision frequency and a trapezoidal electron density 
distribution like that shown in figure 4 were given in reference 2 . A comparison 
of the machine -computed coefficients with those of reference 2 is shown in fig- 
ure 5 as a function of the trapezoid base width. The agreement is excellent, 
but the comparison is unfortunately limited to normal incidence. 
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TABLE I.- COMPUTER RESULTS FOR THE HOMOGENEOUS PLASMA SLAB 


e. 

T 1 

T 2 

R 1 

r 2 

5 rl> 

6 r2> 

s ti> 

5 t2> 

8 r2 ‘ 5 rl> 

5 t2 " s tl> 

deg 





radians 

radians 

radians 

radians 

radians 

radians 




( 

a ) For = 0} z 0 

= 6 Aq - 18 

03 

cm: — = 
o>p 

3 

f 



0 

1.0000 


0 


1.558 


-37.712 










15A15 

.9094 

0.9356 

.4160 

0.3531 

-.0119 

-3.154 

-38.253 

-38.254 

-3.142 

-0.0003 

18.823 

1.000 

1.0000 

.0059 

•0039 

1.558 

4.700 

-32.554 

-32.55^ 

3.142 

.0008 

22.867 

.8735 

.9621 

.4868 

.2729 

-.0119 

3.128 

-33.177 

-33.179 

3.140 

-.0012 

26.161 

1.0000 

1.0000 

.0087 

.0032 

1.554 

-1.585 

-33.854 

-33.852 

-3.139 

.0021 

28.954 

.8216 

.9867 

.5700 

.1628 

-.oi4o 

-3.158 

-34.572 

-3 1 +.575 

-3.144 

-.0028 

31.359 

.9999 

1.0000 

.0145 

.0017 

1.547 

4.694 

-29.074 

-29.069 

3-147 

-.0049 

33.449 

.7471 

.9999 

.6647 

.0106 

-.0160 

3.120 

-29.900 

-29.906 

3.136 

-.0054 

35.261 

.9998 

1.0000 

.0216 

.0018 

1.541 

1.550 

-30.813 

-30.8 o 4 

.0090 

.0090 

36.834 

.6422 

.9803 

.7665 

•1975 

-.0188 

-.0287 

-31.763 

-31.773 

-.0099 

-.0099 

38.170 

.9991 

.9999 

•0455 

.0105 

1.518 

1.538 

-26.550 

-26.529 

.0207 

.0207 

39.287 

.5005 

.8851 

.8658 

.4654 

-.0153 

-.0270 

-27.623 

-27.635 

-.0117 

-.0117 

40.203 

.9943 

.9995 

. 1068 

.0381 

1.453 

1.509 

-28.911 

-28.855 

.0556 

.0556 

40.905 

.3187 

.6487 

.9479 

.7610 

-.0185 

-.0376 

-30.082 

-30.101 

-.0191 

-.0191 

4-1.407 

.9277 

.9858 

•3733 

.1677 

1.179 

1.381 

-25.526 

-25.324 

.2022 

.2022 

41.713 

.0974 

.2175 

.9952 

.9761 

-.0225 

-.0503 

-26.594 

-26.621 

-.0277 

-.0277 




(b) For — 
COp 

= 10" 2 > Z 0 = = 

18 cm; — 

0>p 

3 

~0 



0 

0.7775 


0.0019 


0.0728 


-37.711 










13.415 

.8106 

0.8515 

•3737 

0.3163 

-.0203 

-3.162 

-38.250 

-38.251 

-3.141 

-0.0010 

18.823 

.8624 

.8690 

.0586 

.0390 

.0529 

3.198 

-32.552 

-32.552 

3.145 

.0003 

22.867 

.7683 

.8567 

.4526 

.2398 

-.0207 

3.124 

-33.173 

-33.177 

3.144 

-.0033 

26.161 

.8354 

.8525 

.0815 

.0306 

.0500 

-3.080 

-33.850 

-33.850 

-3.130 

.0007 

28.954 

.7112 

.8535 

.5005 

.1396 

-.0229 

-3.151 

-34.566 

-3^.573 

-3.128 

-.0068 

31.359 

.7952 

.8285 

.1179 

.0147 

.0495 

3-240 

-29.068 

-29.067 

3.191 

.0014 

33.449 

.6340 

.8117 

.5761 

.0101 

-.0251 

3.592 

-29.891 

-29.903 

3.617 

-.0117 

35.261 

.7316 

.7902 

.1777 

.0161 

.0301 

-.0522 

-30.801 

-30.800 

-.0824 

.0011 

36.834 

.5308 

.7527 

.6561 

.1571 

-.0282 

-.0918 

-31.747 

-31.766 

-.0636 

-.01 88 

38.170 

.6240 

.7214 

.2825 

.0790 

.0136 

-.0200 

-26.518 

-26.519 

-.0337 

-.0005 

39.287 

.3968 

.6277 

.7338 

.3528 

-.0272 

-.0805 

-27.586 

-27.612 

-.0533 

-.0256 

40.203 

.4371 

.5779 

• 4677 

.2196 

-.0229 

-.o6o4 

-28.796 

-28.809 

-.0375 

-.0128 

40.905 

.2287 

.4017 

.8042 

.5571 

-.0358 

-.0985 

-29.966 

-30.006 

-.0627 

-.0399 

41.407 

.1679 

.2825 

• 7469 

.5307 

-.0891 

-.1836 

-31.143 

-31.215 

-.0945 

-.0727 

41.713 

.0682 

.1305 

.8704 

.7323 

-.0898 

-.2093 

-25.628 

-25.726 

-.1195 

-.0983 




(c! 

For — 

CUp 

= 10 -1 } z c 

= 6 Aq = 18 cm; -£■ = 
“P 

3 

' 0 



0 

0.0914 


0.0326 


-0.1332 


-37.589 










13.415 

.2758 

0.2800 

.2332 

0.1957 

-.0982 

-3.231 

-38.118 

-38.122 

-3.132 

-o.oo4i 

18.823 

.2526 

.2596 

.2205 

.1508 

-.1393 

-3.258 

-38.683 

-38.692 

-3.118 

-.0088 

22.867 

.2245 

.2557 

.2708 

.1466 

-.1107 

3.076 

-33.005 

-33-020 

3.187 

-.0145 

26.161 

.1986 

.2116 

.2706 

.1086 

-1479 

3.078 

-33.650 

-33.672 

3.226 

-.0214 

28.954 

.1693 

.1855 

.5182 

.0874 

-.1287 

-3.107 

-34.339 

-34.369 

-2.978 

-.0294 

31-359 

.1420 

.1586 

.3330 

.0521 

-.1617 

-2.937 

-35.068 

-35.108 

-2.775 

-.o4oi 

33.449 

.1135 

.1510 

• 3784 

.0288 

- . 1546 

4.160 

-29.557 

-29.609 

4.315 

-.0517 

35.261 

.0875 

.1036 

.4088 

.0482 

-.1845 

-1.008 

-30.359 

-30 . 428 

-.8236 

-.0686 

36.834 

.0630 

.0774 

.4542 

.0964 

-.1930 

-.7193 

-31.190 

-31 . 278 

-.5263 ' 

-.0872 

38.170 

.0426 

.0540 

.4972 

.1526 

-.2232 

-.6653 

-32.013 

-32.127 

- . 4420 

-.1133 

39.287 

.0266 

.0350 

.5430 

.2149 

-.2494 

5.615 

-26.512 

-26.655 

5.865 

-.1430 

40.203 

.0157 

.0212 

.5880 

.2793 

-.2850 

-.7094 

-27.202 

-27.381 

- . 4244 

-.179^ 

40.905 

.0092 

.0128 

.6 261 

.3370 

-.3212 

-.7643 

-27.733 

-27.949 

-.4431 

-.2162 

41.407 

.0059 

.0082 

.6548 

.3823 

-.3533 

-.8176 

-28.090 

-28.338 

-.4643 

-.2482 

41.713 

.0045 

.0061 

.6724 

• 4io8 

-.3756 

-8557 

-28.289 

-28.559 

-.4801 

-.2701 
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-160 
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Figure 4.- Trapezoidal distribution of elec- 
tron density (from ref. 2). 




z o/ x o 


Figure Phase difference of reflected 
components (Ey and Hy) for the homoge- 
neous plasma slab, z q = 

(cup/a)) 2 = 5/9. 


Figure 5*- Reflection and transmission coef- 
ficients as a function of trapezoid base 
width (from ref. 2). f = 1C)9 cps; 
v = 8.241 x 10^ sec" 1 ; 

N e( z l) = 9.5273 x lo9 cm-3 j z 1 = 0.25\). 


Practical Application 

The plasma surrounding a reentering body provides an interesting application 
of the program. The dimensions of many shock layers of current interest are com- 
parable to or smaller than telemetry wavelengths; hence, a W.K.B. solution of the 
wave equation is not valid. 

Typical distributions (ref. 10) of electron density and collision frequency 
are given in figure 6 along a line perpendicular to the body. The N e curve 
which continually increases from the shock to the body corresponds to a com- 
pletely inviscid shock layer. The N e curve which has a maximum between the 

shock and the body is a flow profile which includes a viscous boundary layer. 
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The y = 0 plane defines the H-plane (Ep) and the x = 0 plane corresponds to 
the E-plane (E 2 ) of a slot antenna whose major dimension is along the x-axis. 

The plane-wave transmission coefficients are plotted in figures 7 and 8 for 
the two flow-field approximations. The shapes of the two curves are similar. 

The transmitted perpendicular component of E is maximum at normal incidence 
and transmission gradually decreases with increasing angle of incidence. On the 
other hand, the transmitted parallel components of E gradually increase with 
increasing angle of incidence. At near grazing incidence the parallel transmis- 
sion coefficient drops rapidly to zero. It must he emphasized that figures 7 
and 8 are the result of calculations based on an idealized plasma model; that is, 
transmission and reflection coefficients have been calculated for plane mono- 
chromatic waves incident at various angles on a plane parallel plasma slab in 
which the plasma properties vary only in the direction normal to the bounding 
planes. In lieu of a better method, these calculated transmission coefficients 
can be applied to the free-space radiation pattern of an antenna. Whether the 
resultant radiation pattern is a good approximation of the actual pattern depends 
on many factors. Consequently, some caution should be exercised in using this 
approach . 

The computed results for several angles of incidence for the typical plasma 
layers with and without a boundary layer are given in table II. 



Figure 7.- Transmission coefficients in 
decibels as function of angle of inci- 
dence for the inviscid flow field. 


Figure 8.- Transmission coefficients in 
decibels as function of angle of inci- 
dence for the flow field which includes a 
boundary layer. 
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TABLE IX.- COMPUTES RESULTS FOR THE INHOMOGENEOUS PLASMA SLAB 


e. 


t 2. 

1 - nf. 

1 . if. 

6 rl> 

5 r2> 

& tl> 

5 t2> 

6 r2 " 5 rl> 

8t2 ~ 6 tl 

deg 

db 

db 

db 

db 

radians 

radians 

radians 

radians 

radians 

radians 



(. 

a) Invicid flov 

results; 

f = 244.3 x 106 

cps; z 0 

= 3-36 cm 


o 

-54.8 


-22.1 


-3 106 


1 iiAQ 




10 

-35.0 

- 34.9 

-22.1 

-22.0 

-3.106 

0.0367 

-1.410 

-1.417 

3.143 

-0.0072 

20 

-35. 4 

-34.5 

-22.5 

-21.8 

-3.108 

.0384 

-1.410 

-1.411 

3.146 

-.0004 

30 

-36.1 

-33.8 

-22.7 

-21.4 

-3.110 

.0417 

-1.411 

-1.400 

3.152 

.0114 

45 

-37-8 

-32.2 

-23.6 

-20.5 

-3.116 

.0510 

-1.413 

-1.372 

3.167 

.0411 

6o 

-4o.8 

-29 A 

-25.2 

-18.9 

-3-124 

.0720 

-1.416 

-1.320 

3.196 

.0962 

70 

-44.1 

-26,5 

-26.9 

-17.1 

-3.129 

.1051 

-1 . 4l8 

-1.250 

3.234 

.1682 

85 

-56.0 

-17.6 

-32.9 

-10.8 

-3.138 

.4001 

-1.422 

-■7919 

3-539 

.6239 

89.5 

-76.0 

-18.9 

-42.9 

-7.9 

-3.1+1 

2.199 

-1.423 

.5358 

5-340 

1.959 



(b) Boundary -layer results; 

f = 244.3 x 106 cps; z 0 - 3.36 cm 


0 

-18. 1 


-15.7 


-P Q77 






10 

-18.2 

-18.0 

-13.8 

-13.8 

+- • y / f 

-2.979 

0.1592 

--L . 

-1.340 

-1.342 

3.138 ' 

-0.0026 

20 

-l8.6 

-17.6 

-14.1 

-13.6 

-2.986 

.1666 

-1.345 

-1.330 

3.153 

.0152 

30 

-19.3 

-17.0 

-14.5 

-13.1 

-2.998 

.l8o4 

-1.354 

-1.308 

3.179 

.0411 

45 

-21.1 

-15. 4 

-15.6 

-12.0 

-3.024 

.2196 

-1-373 

-1.249 

3.244 

.1239 

60 

-24.0 

-12.9 

-17.5 

-10.1 

-3.058 

.3052 

-I.399 

-1.133 

3-364 

.2652 

70 

-27.3 

-10.3 

-19.4 

-8.0 

-3.085 

.4323 

-1.418 

-.9766 

3.517 

.4415 

85 

-39.1 

-5.4 

-25.9 

-3.5 

-3.127 

1.249 

-1.450 

-.1782 

4.376 

1.271 

89.5 

-59.1 

-15.6 

-36.1 

-8.7 

-3-1+0 

2.807 

-1.459 

.9538 

5.947 

2.415 


SUMMARY EQUATIONS FOR MACHINE CALCULATIONS 


Accurate solutions of the wave propagation equations for a one-dimensional 
plane-parallel plasma slab can be obtained by numerical integration of the fol- 
lowing differential equations: 

Eq case: 


— + ko 2 (V*r - Ws) = 0 
dz 2 


~ + kQ 2 (Wr + V*s) = 0 
dz 2 
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E 2 case: 


d 2 t 

V 

/dV 

dt _ 


du\ 

w 

/dW 

+ 

dV 

OJ 

IS5 

V 2 + w 2 

\dz 

dz 

dz 

dz/ 

V 2 + W 2 

Idz 

dz 

dz 


V*t - Wu = 0 


d 2 u 

dz 2 


V 

fdW 

V 2 + W 2 

idz 


dt 

+ 

du\ 

+ w _ 

(&Y 

dt 

dW 

dz 

dz 

dzi 

V 2 + w 2 ' 

idz 

dz 

dz 


Wt = 0 


V* = V - sin 2 0 


V = 1 



W = 



At normal incidence only the first set of equations need he solved. The 
starting point for the solution to the equations is z q, and the following bound- 
ary conditions must he met: 


r (Z 0 ) = 1 

s.(z 0 ) = 0 
r’(z 0 ) = 0 
s'(zo) = Eq cos e 
t(z 0 ) = 1 
u(z 0 ) = 0 
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t'(zo) = -kQW(zo)cos 0 


u'(z 0 ) = koV(z 0 )cos 0 


The bo un dary conditions at the surface z = 0 determine reflection coeffi- 
cients, transmission coefficients, and phase factors such that: 




6t2 = Transmitted phase = arc tan 


(-Eg) 

a 2 


- kQZQ cos 0 


A X = r(0) + S ' ^ — 
k q cos 0 


Bi = s(0) - 


r'(0) 

kQ cos 0 


Cl = r(0) - 


s * (o) 

kQ cos 0 


D]_ = s(o) + r ' ^ — 
]£q cos 0 


A 2 = t(o) + 


-W(o)t'(Q) + v(0)u'(0) 
kQ cos 0 W 2 (0) + V 2 (0)j 


B 2 = u(0) 


V(0)t»(0) + W(0)u'(0) 
IkQ cos 0 W 2 (0) + V 2 (0)j 


C 2 = t(0) + j W( 0 )t * ( 0 ) - Y(0)u'(0) 
IkQ cos 0[W 2 (O) + V 2 (0)J 


D 2 = u(0) + 


vCcQfCo) + w(o)u* ( 0 ) 
kQ cos 0 W 2 (d) + V 2 (0) 
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For arbitrary orientation of the electric vector with respect to the plane 
of incidence. 



1 T | 

- (* 5 2 + * 





t R 1 

- (a* 2 ♦ * ?t' k 



Aj = Tj/^cos 20 

+ Tg 2 sin 2 0 

cos 2| 

[ 5 t2 ' 

■ Stl) 

B 3 “ 

T 2 2 ! 

sin 20 sin 2 | 

(^2 " 

5tl) 


A|^ = R-^cos 2 ^ 

+ R 2 2 sin 20 

cos 2| 

k 6 r2 ‘ 

' s rl) 

% = 

R 2 ^sin 20 sin 2 | 

( 

, s r2 " 

s rl) 



If the incident wave is elliptically polarized, the polarization phase factor £ 
is added to and &tl- The solutions to the wave equation remain unchanged. 


CONCLUSIONS 


The equations and the associated boundary conditions which describe the 
interaction of electromagnetic waves with a one-dimensional inhomogeneous plasma 
slab were developed in a form suitable for machine programing. The program solves 
for the reflection coefficients, transmission coefficients, and phase factors for 
the components of the electric vector perpendicular to and parallel to the plane 
of incidence. From these results, a wave for which the electric vector is of 
arbitrary orientation to the plane of incidence and of arbitrary polarization can 
be completely defined. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., April 23 , 1963. 
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APPENDIX 


IMPLICATIONS OF BOUNDARY CONDITIONS 


The purpose of this appendix is to show: 

(1) That the behavior of the wave in the plasma as a function of x is of 

ikQXsin0 

the form e 

(2) That the exit angle of the wave is equal to the angle of incidence 

(3) That the solutions of the wave equations are independent of the 
y-coordinate if the incident waves are independent of y. 

The general form of equation (6) implies solutions which are functionally 
dependent upon all the coordinates x, y, and zj that is, 


Ey = Ey(x,y,z) 


(Al) 


If by arbitrary choice of coordinates, the incident and reflected waves in free 
space are independent of y, then the solution of the wave equation is 


Ey(x,0,z) 


e ik 0x x 




(A2) 


The term is dropped, so that the waves are required to propagate 

in the positive x-direction. The coefficients M and N are interpreted as 
incident and reflected wave amplitudes, of such character that: 


Eiy(x,0,z) = Me ik OxX+iko zZ 

= E iOy e ik O xsin0+ik O zcose (A3) 

EpyUjOjz) = Ne ik 0x x - ik 0z z 

= Eroye 1 ^ 81110 " 11 ^ 20080 (a4) 

where, from figure 1, kQ X = Isq sin 9, kQ Z = kQ cos 0, and 0 is the angle of 
incidence. 
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The boundary conditions require that the tangential components of E he 
continuous across the surface z = 0. The boundary condition for the perpendic- 
ular component is: 


— » 
Ury 


% u y 


= u„ 


fs 


%Uy 


(A5) 


After the vector operations are performed, the boundary condition becomes 

^yf s = -^yp (A6) 

The value of Eyf s is the sum of equations (A3) and (a 4) . Therefore, at 
z = 0, 

Ei0ye lk 0* sln9 + Eroye 1 ^ 51110 - Ey(x,y,0) (AT) 

At x = y = 0, equation (A7) becomes 

■^iOy + ^rOy = £y(0j0,0) (A7a) 


The plasma properties at a point on the surface z = 0 are indistinguish- 
able from those at any other point. Hence, equation (A7a) must be satisfied at 
any point on the surface z = 0. At another point (0,y2)j equation (A7) is 


E i0y + -^rOy ~ -%"( 0 ^y 2 ^ 0 ) (A7b) 

Equations (A7a) and (A7h) are identical. Therefore, the electric vector in the 
plasma is independent of the coordinate y, and vave propagation must occur in 
the xz -plane. 

At another point (x 2 , 0 ) equation (A7) becomes 


(EiOy + E r0 = Ey(x 2 ,0,o) (A7c) 
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The presence of the sinusoidal factor e lk 0 x 2 sin implies that the energy trans- 
mitted into the plasma varies as a function of x. This implication is physically 
unreasonable for incident plane waves , and the right-hand side of equation (A7c) 
must be 


%(x 2 ,0,0) = % (0,0)e lk 0 X 2 8in9 
Thus , the right-hand side of equation (A7) has the form 

Ey = E y ( Z )e ik O Xsin9 


S F v (z)e ik O XSln0 


(A8) 


(A9) 


The solution to the wave equation in free space for z ^ zq is 


E ty _ E tOy e 


i 


k tOx x+lk tOz z 


(A10) 


-iiu.Q Z z 

An e term does not appear because the wave is outgoing. If 0-^ is the 

angle between the z-axis and the direction of propagation, then 


E ty ~ E tOy e 


ikQXSinGt+ikQZcose-t 


(All) 


One of the necessary boundary conditions requires that equations (A9) and (All) 
must be equal at z = zq • Therefore , 


F y (z o )e ik O xsin0 = Et0ye ik 0 xsin 9t +ik 0 z 0 cos ®t 


(A12) 


At x = 0, equation (A12) becomes 


F y( z o) - E t 0 y eik ° Z o COSe t 


(A12a) 


Since the plasma properties are indistinguishable over the surface, equa- 
tion (A12a) must be identically satisfied for every point (x,y) at z = zq. 
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Therefore, equation (A12a) completely specifies Fy^ZQ^. Substituting equa- 
tion (A12a) into equation (A12) leads to the identity 


^ikQXsine = ^ikQXsinQ-t 


(A13) 


Hence 


e = e t 


The continuity requirements for the component Hy, lead to the following equation: 

^H i0y + H r0y j e ik O XSin0 = H y (x,y,0) (Al4) 

This equation is of the same form as equation (A7)> and arguments similar to 
those applied in the analysis of Ey show that 


Sy(x,z) = Gy(z)e lk O xsine 


(A15) 


It is of interest to note that 
solutions , Ex , Ey, E^j Hx^ 


Hy is also independent of y; hence all the 
Hy, and H z are independent of y. 
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